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Es HISTORY AND PURPOSE 


A. SATELLITES 


The space age kicked off in 1957 with the Russian 
launch of Sputnik. Sputnik-1l was immediately followed one 
month later with Sputnik II, which carried a dog named 
Laika. The United States countered with the launch of 
Explorer 1 in 1958, which went on to discover the Van Allen 
radiation belts. These feats led to the launch of Echo 1, an 
aluminum-coated plastic sphere that passively reflected 
voice and picture signals. Next was the first true 


communications satellite, Telstar. Telstar actively relayed 





communications between locations that were at great 
distances from one another [1]. Over the years, the use of 
satellites has evolved, and now they function in many 
different roles. Satellites are used for navigation, as 
relay stations for radio and television, communications, 


weather, observations of the Earth and space, and for wars. 





The use of these satellites will not diminish; rather, as 





time progresses, there will be new inventive things that can 


be done with a satellite. Additionally, as resources become 





scarcer, the satellite will need to get smaller and cheaper. 


The Galileo, Cassini and Hubble spacecrafts had life 


cycle costs over one billion dollars. The Clementine space 





craft cost $85 million and the Lunar Prospector cost $63 


million. After the Cold War ended, the exclusivity of space 





to a select few countries started to change. There became a 
commercialization of space, from the previously 


military/government specific one. There was an increase in 


the availability for launch services and a dedicated launch 


vehicle was no longer required. It became cost effective to 


get as many different companies, or entities, on 


one launch 





satellites, 





vehicle as possible to offset the costs. Large 


such as in Figure 1, can cost from $5,000 to $20,000 per 


kilogram to launch, and from $30,000 to $160,000 per 





kilogram to build [1]. One driving force in the design cost 





for satellites are the use of space-qualified components. 





Figure l. JCSAT-12 being built. From [2] 


Almost all companies, including the U.S. 





government, 


require the components of a space craft or satellite to 


contain space-qualified components. These components must 


meet a military or NASA (National Aeronautics 


Z 


and Space 


Administration) standard for the material, fabrication, 


storage and transport. Thes requirements significantly 








increase the price of a satellite. Additionally, anywhere 








from two to four models must be made; one for testing, one 





for “fight, a possible flight spare and a_e possible 
engineering model [1]. Costs need to be cut wherever 


possible. Decreasing the siz requires more advanced 








technology that is smaller than its predecessors. Satellites 


need to be lighter and more efficient. The use of commercial 





off-the-shelf (COTS) products dramatically reduces these 
costs. Multiple spares can be purchased for a fraction of 
the price of designing and fabricating required circuitry. 
These steps have led to an increase in the use of satellites 


by universities for education and research [1]. 


Satellites are classified by weight into size 
categories. Cube Satellites, Figure 2, which are known as 


CubeSats, can be nano-satellites (weighing from 1-10 kgs) or 








pico-satellites (weighing less than 1 kg) and are generally 











10 centimeters cubed. They are orders of magnitude smaller 


than their full satellite cousins. The idea of a CubeSat is 








not only important to the space industry, but to education 





as well. CubeSats are becoming important and useful to high 


schools, colleges, and universities. Students are designing, 





building and launching these small objects through their own 
programs. The knowledge they gain in these programs prepares 
them for joining the industry’s satellite community. Their 


small size allows CubeSats to be placed as a_ secondary 





payload on an already scheduled launch vehicle, vice having 
the dedicated launch required for larger satellites. For the 


industry, a CubeSat can be launched on a budget of only 





$100,000 and the launch can take place years before a full- 
3 


sized satellite could [3]. For student programs, there is a 
discount to stimulate research and development for new 
innovations, where a CubeSat can be launched for as little 


as $40,000 [4]. This price allows new and radical 





technologies to be tested in space that would not normally 
take place due to the large monetary investment that the 
industry has on a given launch. Industry needs systems that 
are guaranteed to work and do not have the money to spend on 
a “what if.” These CubeSats allow for a much lower cost of 
space testing for these advancements. Launch opportunities 


still remain an issue. Two schools, California Polytechnic 





State University (Cal Poly) in San Luis Obispo, CA, and the 
Naval Postgraduate School (NPS) in Monterey, CA, have been 


instrumental in assisting in this area. 











Figure 2. Size of a Cube Satellite. From [5] 


B. STUDENT INVOLVEMENT 


Cal Poly has standardized the industry by creating the 
Poly Picosatellite Orbital Deployer, or, as it is commonly 


known, the P-POD. The P-POD (Figure 3) can hold up to three 


units (3-Us) of CubeSats. A single unit of CubeSat for the 
P-POD is defined as a 10 cm x 10 cm x 10 cm satellite. As 


such, a two unit (2-U) satellite would be 10 cm x 10 cm x 20 





cm. Each P-POD contains and shelters the _ satellite(s) 
throughout the launch process and once orbital altitude is 
reached, deploys the cargo. The P-POD system is self 
contained and standardized. CubeSats going into the P-PODs 
have maximum weight limits and center of gravity 


constraints. Once CubeSats are built, they are sent to Cal 





Poly, where they are incorporated into the P-POD and sent to 
the launch site for transportation into space via a 


previously arranged launch vehicle [6]. 





Figure 3. Cal Poly Mk II P-POD. From [7] 


A common problem in satellite launches is missed launch 
opportunities. Parts not being delivered on time, lack of 
adequate access to clean rooms, failed inspections/tests, 


and countless other reasons also can cause the delivery date 





of satellites to be missed. Current rocket launch procedure, 
the method of CubeSat delivery to space, requires a hard and 


fast deadline for delivery of all cargo. This deadline is 


essential to calculations of fuel and launch criteria for 





the rocket. If a CubeSat/P-POD fails to be delivered on 








time, weight equivalence is used and it is no longer 





possible to get onto that flight. The flight opportunity is 








missed and must be rescheduled. Flight opportunities are 


rare and that is why the P-POD is crucial. This allows three 








CubeSats to be launched instead of one. 


Cc. NPSCUL 


Currently, there are about 40 to 50 university programs 


developing CubeSats [4]. All of these programs would like to 





test out their systems in space but, even with Cal Poly’s 





three-unit launcher, there are still not enough 





opportunities. As such, NPS has stepped in to assist with 


the creation of the NPS CubeSat Launcher, or NPSCuL as seen 





in Figure 4. NPSCuL, in short, is a unit that contains eight 


P-PODs, increasing a rocket’s delivery potential from 3-Us 





of CubeSats to 24-Us. Therefore, when one launch opportunity 
comes along, multiple educational programs can get on board. 


This can also lower the individual cost per launch. 





Figure 4. Image of NPSCuL. 


NPSCuL has already been designed, built and 


qualification tested to ULA ABC launch environment [8]. 





NPSCuL was designed specifically for ESPA (EELV or Evolved 
Expendable Launch Vehicle secondary payload adaptor). The 
ESPA ring (Figure 5), allows for six payload attachments. 
NPSCuL can be connected to ESPA compatible fixtures on the 
Atlas, Athena, Delta, Falcon 9, Falcon le, Minotaur IV, 
Taurus and European Space Agency’s VEGA rockets [9]. On the 
Atlas V, one possible union of the launcher and rocket is 
through the United Launch Alliance (ULA) Aft Bulkhead 
Carrier (ABC) [10] (Figure 6). The NPSCuL structure consists 
of a 5-sided box, 8 P-Pods, a sequencer, and fasteners. All 
specifications for the current design of NPSCuL, as tested, 


are shown in Table l. 





Primary Payload / Satellite 





Figure 5. ESPA ring on the EELV. From [13] 








Existing 
Mounting 
Brackets 


Centaur Aft 
Bulkhead 


Aft Bulkhead 
Carrier (ABC) 


Secondary 
Payload (SP) 


Figure 6. ULA ABC. From [14] 


Table 1. Current NPSCuL structure data. From [11] and [12] 









































NPSCuL lbs 

8 P-Pods 116.84 
Sequencer plus fasteners 8.80 
Box plus fasteners 42.15 
Harness 2.90 
Total weight of NPSCuL 170.63 
Max allowable weight of structure 177.00 
Current usable weight 6.37 





Even though there is a current design for the NPSCuL, it 


could use some improvement. While technology is getting 





smaller and better, it also allows satellite designers to put 














more equipment into the CubeSat structure. With the 
advancement in technology making parts smaller, developers 
want to fit more into the CubeSat. As such, this increases 


the demand of the previously allotted 1.33-kg weight 











allowance [6]. By making the structure lighter, this will 








allow the freed weight to be applied back into the CubeSats. 








The design constraints are shown in Table 2. The weight 
constraint for launch still remains with a maximum weight of 
NPSCuL, including the CubeSats, not exceeding 177 lbs or 
being lighter than 169 lbs [11]. 





Table 2. Design constraints for NPSCuL system. 
From [11] and [15] 













































































Max weight of structure 177.00}1lbs 
Minimum weight of structure 169.00|;1lbs 
Weight of a side of the box 7.44/1lbs 
Expected g-forces in X direction [12] 5|g's 
Expected g-forces in Y direction [12] 5ig's 
Expected g-forces in Z direction [12] T\g's 
Density of Al 7075-T7351 [11] O4.101.| tbs /aint3 
Modulus of elasticity [11] 10.3E6|/psi 
Yield Stress [11] 57) ksi 
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II. THE DYNAMIC STRUCTURE 


A. LAUNCH EFFECTS 





NPSCuL must get to space to be effective. The current 





method is to go via rocket and launching on a rocket is not 


a smooth ride. There are strong forces, up to seven times 





the gravitational constant, and random vibrations. The 


natural frequencies of the structure are of the utmost 











interest. At these natural frequencies, the CubeSat launcher 


can see a large amplification of forces, which in turn will 











be transmitted to the CubeSats. Even though the launcher 


structure can easily withstand these forces, the safety of 





the CubeSats inside the launcher becomes a question. The 


satellites, though built tough for the space environment, 





still have sensitive sensors that can be damaged. The 


current design for NPSCuL has already been analyzed using a 





Finite Element Model (FEM) in NX I-DEAS 6; a modeling 





software. The first six modes of the structure are selected 


for reference in this research. 


B. MODAL SHAPES 





The modal shapes refer to the movement of the structure 


at the system’s specific natural frequencies. Of NPSCuL’s 





lowest six modes, this analysis will only show modes 1, 2, 3 





and 6. Mode 1 is the lowest natural frequency, and modes 2, 
3, and 6 have significant center of gravity movement and 


structural displacements. Figures 7 through 10 show the 














analytical display of the analysis of these frequencies [17]. 
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Figure 7. Mode 1 analyses. From [17] 





Figure 8. Mode 2 analyses. From [17] 
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Figure 9. Mode 3 analyses. From [17] 














Figure 10. Mode 6 analyses. 


From [17] 
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Mode 1 is considered “corner flexing” [17] and is the 





lowest natural frequency for the structure. The requirements 
for the structure is to have a first natural frequency no 
lower than 35 Hz [15]. The first natural frequency was 
analytically evaluated at 41.39 Hz. Figure 7 shows the 
rendering of movement for mode 1. Mode 2 and mode 3, 53.75 
and 54.73 Hz, respectively, are modes that are a “rocking 


motion” in the X-Y plane, as seen in Figure 8 and the Y-Z 





plane, as seen in Figure 9. Mode 6 is a “po-go” motion of 
the entire structure up and down onto the adaptor ring. This 


is shown in Figure 10 [17]. 


Cc. VIBRATIONS TESTING 


The next step is to take the results determined by 
analytical methods and verify them with a vibrations test 
(Table 3). The engineering model of NPSCuL was loaded with 


mass models for the P-PODS and, with the assistance of tri- 





axial accelerometers, the actual natural frequencies wer 





found. Table 4 shows the results of this test [10] and Table 5 
shows this data compared to what was expected from the 
analytical results. Figures 11-14 show the structural movement 


of the first and primary modes (modes 1, 2, 3 and 6). 
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Table 3. Comparison of I-DEAS modeling and test results. 


Analytical Tested Percent 


[a ea.ce [03.93 [ 0.89 
[6 [95.06 | 100.5 [ 5.41 





DEFORMATION: 6-6: COMBINED_100403_sDOF_sSHAPE/42 95339 
MODE: 6 FREQ: 42.9589 Damp: 0.5224496 
Acceleration - MAG MINH: 0.006400 max: 4.37&-00 

FRAME OF REF: PART 





Figure 11. First natural frequency, 42.96 Hz, structural 
deflection. From [10] 
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DEFORMATION: 4-4: COME INED_100408_s00F_ phe © 83442 
MODE: 4 FREQ: 57. 83441 camp: 1.987136 
Acceleration - MAG MIN: 0,00E+00 max: 2. 52e601 

FRAME OF REF: PART 





Figure 12. “Rocking” natural frequency, 57.83 Hz, structural 
deflection. From [10] 


DEFORMATION: 7-7: COMETNED_100408_s0OF_sHaPE /59.32238 


MODE: 7 reed: $9.32238 OOMP: 1.93505 
Acceleration - MAG MIN: 0. 008.00 max: 1.29.01 
FRAME OF REF: PART 





Figure 13. “Rocking” natural frequency, 59.32 Hz, structural 
deflection. From [10] 
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DEFORMATION: 2-2: COMBINED_100408_s00F_sHAPE /100.4971 
MODE: 2 FREQ: 100.497 Dame: 0.2534246 
acceleration - MQG MIN: 0.008400 max: 7.15.00 
FRAME OF REF: FART 





Figure 14. “Po-go” natural frequency, 100.50 Hz, structural 
deflection. From [10] 


D. OBJECTIVE 


The goal is to reduce the mass of the current 
structure. The NPSCuL is designed to create a way to launch 
more CubeSats per launch opportunity. The design shape was 
determined by fitting the structure into the given volume 
and mass constraints. The following analysis is the first 
step in optimizing the mass of the structure to see where 
weight can be removed without unduly reducing its strength. 
Once a reduced weight is achieved, further analysis to 


include dynamic response can be made with the new design. 


There are two possible directions for weight removal: 
modify the current structure, or completely redesign it. The 
first option of modifying the structure is the cheapest and 
fastest way as the time, materials, and manpower do not have 
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to be spent to start over. An analysis needs to be done on 





the weight of the system to see if weight can be removed. 


This will be the initial iteration at looking at 








modifications to the structure. Once weight can be removed, 








it can either be integrated back into the CubeSats’ weight 
allowance or remain off. Additional alternatives can include 
external stiffeners, a damping system, or modifications to 
the adapter ring. All of these options eventually go toward 
reducing the total mass of the launcher or, in the case of 
the alternatives, will reduce the forces seen in the 


CubeSats during launch. 


The second option would be to change the structure from 





its current box configuration and to start with an entirely 
new design. Using a different shape or reducing the number 
of CubeSats the structure can hold are all possible 


alternatives, but are beyond the scope of this analysis. 
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IIIT. SIMPLIFICATION OF THE STRUCTURE 


A. SINGLE WALL 


The first simplification needed is to take the four- 





walled box and look at one single wall. As seen in the 
Figure 15, a single wall has two P-PODS attached onto the 
plate on the right-hand side. This area is not going to be 
modified since nearly a quarter of the weight of the 


structure is concentrated in this region and there is a 





large presence of fasteners. The dots represent the 





locations of the fasteners for the P-PODS, and the area to 


be modified is the 6.365-inch by 12.0-inch section on the 





left. For this analysis, the wall depth remains at a uniform 


0.25 inches throughout the wall. 






12.0 Inches 
S2y2u] 0 €T 





6.365 Inches 





Figure 15. Two-dimensional view of side wall, including the 
modifiable area (MA). 
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B. COMPARISON TECHNIQUE 


Modifying the wall volume of the current structure, 





while maintaining all remaining dimensions, allows an 





analytical comparison to the original test data. The 





original structure, of a solid Al 7075 wall plate with 
dimensions of 18 inches by 13 inches by 0.25 inches, is used 
as the baseline when comparing the total mass of a full 
wall. When comparing solutions for the modifiable area (MA), 
the baseline is an Al 7075 wall plate with dimensions of 
6.365 inches by 12 inches by 0.25 inches. Analyzing the 
desired alternatives would show an increase or decrease in 


mass and/or an increase or decrease in stiffness. The 





following analysis determines the baseline of the MA by 
using Al 7075. 


Cc. FLAT PLATE ANALYSIS 


To analyze the stress and deflection seen in the MA, 





the gravitational forces must be turned into a distributed 





load across the top of the wall. Once the equivalent force 








is calculated, the stress and deflection can be determined 
and from there the stiffness can be calculated. The maximum 


weight of the structure is 177 lbs. During launch, the 





structure can see up to 5 gs of force in the X-direction, 5 
gs of force in the Y-direction and 7 gs of force in the Z- 
direction. Using the static-equivalent to these dynamic 


forces, the following calculations show why 1200 lbs is 








used. The 6.365 inch by 12 inch MA is approximately one- 


third of the total volume of one side wall. The factor of 





safety used for this calculation is two. 
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1 gravitational unit = 9.8067 N/kg = 1.0000 lb¢/lby 





Combined load = V(5°+5°+77) = 9.95 g’s = 9.95 lb¢/1bm 
9.95. Ubefiby * TTT Thy AOL Let ee 
T7615 Tbe * (1/3) = 587.05: 1b: 


587.05 lbs * 2 = 1200 lbe 


Now, the full unmodified wall can be analyzed. The 





current structure is made from Aluminum 7075-T7351 and a 
full plate analysis yields the baseline data for comparison 


to volume removal in the modifiable area. The material 





properties are found in Table 4. 


Table 4. Aluminum 7075-T7351. From [18] 





Modulus of Elasticity (E) |10.3 x 10° | psi 








Density 0.101 lbm/ in? 
57 ksi 


Yield stress (oy) 





68 kea 


Ultimate stress (0y) 

















The mass of a single full wall is equal to 





p*t*l*h = 0.101 ibe/ine -* 0.256 in * 164m * 13° Aine =-5..909- Iby 
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The mass of the modifiable section of a single wall equals 
pFt*l*h = 0.101 Ibs/in® * 0.25 in * 6.365 in * 12 in 
= 1.929 lb» 
Maximum stress is [16] 


o= PVA = 1200 Tbe 7 (0,25 am * 6.365 an) = 754,12 psi 


The stiffness in the area to be modified [16] 
k = [AE]/l = [(0.25 in * 6.365 in) * 10.3x10° psi] / 12 in = 
1366x10° lb-/in 





These values will be essential in further calculations and 
show that the modifiable area is approximately one-third of 


the total mass of one side of a wall. It is also important 





to keep in mind that one of the four side walls contains the 





sequencer and will not be modified. This wall is shown in 


Figure 16. 


12.0 Inches 





6.365 Inches 


Figure 16. Two-dimensional view of side wall with sequencer 
bolt holes. 
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IV. MASS REDUCTION OF THE ALUMINUM WALL 


A. REMOVAL OF WALL VOLUME 





Due to the presence of corner flexing in the current 





structure during modal movement [10] and [17], simply 


removing the 6.365" by 12" section is unwise. The complete 





removal would further undermine the corner regions, adding 
to the flexing already seen. There need to be diagonal 
stiffeners in this region. The stiffness of this modifiable 


area started as 1366x10° lb-/in; however, onc th center 








section is removed, it becomes approximately 





oS TRET/ LS 110,25 Se PD ay See? per} 7 12: anc = 
214.6x10° lb-/in 








This reduces the original stiffness by nearly 85%. Table 5 








shows the comparison of a full wall of A1-7075 compared to 





removing the entire modifiable area. Figure 17 shows this 


data graphically. 








Table 5. Comparison of removing all of the modifiable area. 
Percent of Percent of 
Mass mass Stiffness stiffness 
Full Wall 5.909 100 1366 100 
Removal of MA 3.980 67.4 214.5 Bag 7 
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Figure 17. Visual display of full wall versus removing the 
MA. 





This drop in stiffness is dramatic. The verification to 





ensure that this stiffness drop does not undermine the 
system is to check for buckling of the thinnest members. 
Checking for buckling along with verifying that the yield 


stress of the material is not violated and the deflections 





are minimal determines the suitability of the alternative. 





The sides of the MA are 0.5 inches wide; the wall has a 





thickness of 0.25 inches, and the member seeing the largest 
force is 12 inches long. Using the buckling criteria [19] 
for a beam, the critical load that can be carried before 


failure is calculated. This value can then be compared to 





the maximum force seen in the member. The results of this 
analysis, as shown in Table 6, explain that the removal of 


the entire MA results in failure. 





Per = [n?*E*I]/L? 
Where I = [w*t?]/12 
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Table 6. Buckling failure of full MA removal. 


EF Pcr 
MA (lbs) (lbs) | Status 





[0.5] 600] 459.6 


B. TRUSS ANALYSIS 


The possible solution to optimizing the mass is for 
material from the MA to be removed while keeping portions as 


support members. To do this, a Finite Element Methods (FEM) 








technique is utilized. FEM is a way to take a complicated 





problem and simplify it by putting it into a computer 
program to be analyzed. Three-dimensional problems can be 
Simplified as two-dimensional ones with the use of beam 
elements. Using the basic MATLAB code found in [20] for a 


truss system, the modifiable area was turned into a truss 








problem. FEAPLYC2.m, FEASMBL1.m, FEELDOF.m and FETRUSS2.m 





were all used as the basic code and then a function file was 
created to run these programs with various truss designs. 
Only the code in the function file was modified for each 
variation in the design. The analysis was conducted on eight 


different truss designs, as shown in Figures 18-20. 
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i’ ii - 


% “| “ 606 Ib 600 Ib 
Cc 
Figure 18. Truss variations A, B and C. 


il ‘id ‘fi 


“| “ “| “| “| 
Figure 19. Truss variations D, E, and F. 
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r r" 
600 Ib 606 Ib 600 Ib 600 Ib 


| 


G H 
——»x ——ik 
Figure 20. Truss variations G and H. 
Cc. CALCULATIONS 


Not only are all eight configurations analyzed, but 
additionally each configuration is evaluated using three 
different truss beam widths; 1 inch, 0.5 inches and 0.25 


inches. For all cases, each beam remained at 0.25 inches of 





thickness and the beam width refers to the width of the 
face. An example of how the three different beam widths 
would look on configuration A is shown in Figure 21. The 
function codes for each variation can be found in Appendix 
A. Each file is designated first by the letter assigned to 


the configuration as shown in Figures 18-20, and secondly by 





the size of the beam width. FEMA1 denotes the FEM analysis 
of a l-inch beam width, FEMA5 denotes the 0.5-inch beam 





width and FEMA25 denotes the 0.25-inch beam width. 
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oy A hy 











600 Ib 600 Ib 600 Ib 600 Ib 
x = i 
1” width 0.5” width 0.25” width 
Figure 21. Example of what different beam widths look like in 
the MA. 


The output created by MATLAB contains the maximum 
stress in each beam element and also each node's deflections 
in the X and Y directions. Elements and nodes are numbered 
in the clockwise direction. Each node has two degrees of 
freedom, the X and Y directions. As such, the first node 
would have the designation of the X direction’s deflection 
denoted as one, its Y direction’s deflection denoted as two 


and node two would have X and Y directions listed as three 





and four, respectively. These results for all 24 
permutations are found in Appendix B. From the _ stress, 
deflection, and volume of removed material, each variant’s 
percent of total wall mass and stiffness can be calculated. 
The following data is used to create the output found in 


Table 7. 


Cross-sectional area of 1" width beam = 0.25 in?’ 


Cross-sectional area of 0.5" width beam = 0.125 in?’ 





Cross-sectional area of 0.25" width beam = 0.0625 in?’ 


Volume of Modifiable Area (VMA) = 0.25 * 6.365 * 12 
= 19.095 in? 
Volume of Voids (VV) = Volume of material no longer present 
in MA 


Table 7 shows the calculations for the stiffness and mass 


loss for each situation. The percent decrease in mass is 





plotted against the percent of stiffness from the original 


solid wall baseline. Figures 22-24 show the results for beam 





width thicknesses of 0.25, 0.5 and 1 inch, respectively. The 





most advantageous system would have the largest percent of 
mass reduction while having the largest percent of 
stiffness. As an example, in Figure 22, configuration A has 
the highest reduction in mass but the lowest amount of 


stiffness. Conversely, configuration C has the largest 





stiffness but is only average in mass reduction. 
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Percent reduction in mass and stiffness calculations. 


Table 7. 


Mass freed 
whole 
structure (Ibs) 


ran 


B) &| B| B 

MGM Mi 

eae oe 
81g 


g = 5 
Ss ale 
5 a Ae 
—— 
CTachchcich 


sid 
sige 
SEEEL 
oe 


ra. | as] one [.006e [ace | azs| 00.5] 1080 Pama | marl asa 


0.25 | -0.0092 | 7.9936 | 7993.6 |0.0625 | 499.6) 54304 
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Figure 22. 


% Reduction of Mass 


Figure 23. 








0.25" Beam Width 
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% of Stiffness 


Reduction of stiffness versus mass for 0.25-inch 
beam width. 
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Reduction of stiffness versus mass for 0.5-inch 
beam width. 
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Figure 24. Reduction of stiffness versus mass for l-inch beam 
width. 


A definitive optimal configuration is not’ easily 
discerned. As mass is the most important parameter for this 


analysis, looking at each configuration’s contribution 





toward total weight removal for the full structure is 
assessed. The computations in Table 7 are used to achieve 
Figures 25-27. The best four choices for all three beam 
width sizes are configurations A, B, C and F. The next step 
is to ensure that none of the configurations buckle in all 
three width sizes. This would eliminate a configuration from 
the analysis. Table 8 shows the calcuations and verification 


to whether buckling occurs. As can be seen, there are three 





cases where buckling would take place. Threfore, 
configuration A can be used only if l-inch beam widths are 
used, and configuration F can be used if 0.5- or 1-inch 


widths are used. 


ae 


Freed mass 0.25" 





Figure 25. Total freed mass from all truss types at 0.25 
inches. 


Freed mass 0.5" 





Figure 26. Total freed mass from all truss types at 0.5 
inches. 


oo 


Freed mass 1" 





Figure 27. Total freed mass from all truss types at 1.0 
inches. 
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Buckling failure for each configuration. 


F=s*A 


(Ibs) Pcr (Ibs) | Status 


Safe 


fe 
3676.83|Safe 
1838.41|Safe_ | 


229.80 


3676.83 


0 
362.38 
8 


1838.41|Safe_ | 


3676.83 


8 
600.00 
0} 1838.41|Safe 


362.3 


0 1 


Now that the best configurations based on weight have 


been established, 
four designs. Figures 
enlarged. 
but the lowest stiffness. 


stiffness but 





Configuration C has the highest 


is tied with configuration B for the lowest 
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it is worthwhile to compare the remaining 
28-30 show the four configurations 


Configuration A has the highest reduction of mass 





mass reduction. That leaves configurations F as a good 


tradeoff between mass and stiffness. 
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Figure 28. Comparison of four favored solutions for 0.25-inch 
beam width. 


0.5" Beam Width 
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Figure 29. Comparison of four favored solutions for 0.5-inch 
beam width. 


36 


1" Beam Width 


22.50 
. 4 


22.00 
21.50 
21.00 
20.50 
20.00 
19.50 
19.00 a 
18.50 + T 1 1 1 1 
0.00 10.00 20.00 30.00 40.00 50.00 























Ww 
2) 
o 
= 
— 
° 
i 
2 
> 
VY 
Ss 
mo] 
i) 
[a 
xs 











% of Stiffness 





Figure 30. Comparison of four favored solutions for l-inch 
beam width. 
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Vv. CHANGE OF MATERIAL 


A. CURRENT MATERIAL 


The current structure is made from space-grade 


Aluminum, Al 7075. This metal is used due to its light 





weight and high resistance to stress corrosion cracking. The 


material is finished with a chemical conversion coating for 








corrosion protection, and to lower the electrical resistance 





es A different approach to lightening the structure is 





to look at changing the metal. An analysis is to be 


conducted looking at not only lighter metals but also 








heavier ones as compared to Al 7075. In order to properly 








analyze a heavier metal, the use of the MA will be required. 





Lighter materials could incorporate the use of the MA but do 
not have to. Firstly, there are six pounds of free weight 
currently available for use in the NPSCuL configuration when 
using the space-grade Aluminum. If a heavier metal is used 
that exceeds this weight constraint, voids will have to be 


made to maintain the system within parameters. From Chapter 





IV, it is seen that for Aluminum, the “F” configuration is 
the most efficient. However, with a change in material, that 


configuration may change. 


B. OTHER MATERIALS 


Using a list of general, known and available metals 
[19], a simultaneous assessment can be made of different 
metals with different beam widths. Appendix C shows the data 
used for the material analyzed. Table 9 shows the range of 


variables input into a computing program called MINITAB. 
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Table 9. Maximum and minimum metallic values for modulus, 
density and beam width. 


























Minimum Maximum 
Modulus of Elasticity (*10° psi) 6.5 (Magnesium) 58 (Tungsten) 
Density (lb/in’) 0.0639 (Magnesium) |0.697 (Gold) 
Beam width (in) O25" 1.0" 




















Cc. MINITAB 





MINITAB is a systems ngineering analysis tool that 


takes a given set of inputs and returns desired outputs. For 








this analysis, the inputs are to reduce the area of material 





in the modifiable section, change the material to something 
other than Aluminum, or a combination of the two. The 


outputs of the MINITAB program are the new mass of the 








changed wall and the corresponding stiffness. These are the 





same outputs seen from MATLAB. However, an improvement of 


MINITAB over MATLAB is its ability to analyze three or more 





parameters simultaneously. In essence, a three-dimensional 


graph can be analyzed by holding one of the variables and 





looking at how the other two variables respond to the held 


value. This is the Response Surface Model (RSM). In the RSM, 





the feasibility solution space is identified. The structure 





is not to exceed 177 lbs. The output graphs will show 





regions of acceptable material Modulus of Elasticity and 








density. From here, a material only needs to be selected 





that fits these criteria to satisfy the constraints. 


Using the information from Table 9, a list of 20 random 
combinations is created. This data is then input into MATLAB 
using configuration F program code. The maximum stress and 
displacements are found and in conjunction with the 


calculated mass, the input is ready for MINITAB. This data 
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is shown in Appendix D. 


Once the mass and stiffness data is 


the program generates the region of acceptance in 





entered, 
the solution. Figure 31 


entered into MINITAB. 


Create Response Surface Design - Factors ee 


Levels Define 
Cube points 





Factor Name Low High 
A M_Elasticity 6.5 58 
B Density .0639 .697 
Cc Truss Width 2 1 


Figure 31. 


The outputs of looking at the desired 0.25-, 


inch beam width 


respectively. In Figure 





inches, 


section in the upper left-hand corner. 





are approximately where the metals Titanium 





elements 


the acceptability region 


rie 


ld be located. 


lb/in*®), Aluminum (densi 
(density of 0.0639) woul 
lowest density material 


shows the constraint of the data 








coc_| 





Parameters used for MINITAB input. 


Q0.5- and 1- 
32-34, 
0.25 


are shown in Figures 





32, where the beam width is 


is located in the white 


Flagged on the graph 
(density of 0.163 
y of 0.101 lb/in*3) and Magnesium 





These are the three 





thickness of beams, 


constraints. Figure 33 h 


the materials still do not satisfy the constraints. 


34, 


s looked at and, for 0.25-inch 


do not satisfy the mass and or stiffness 


as a larger feasibility region, but 


In Figure 


the feasibility region shows that the material needs to 


have a density of less than 0.108 lbs/in® and a Modulus of 


41 


Elasticity between 9.25 x 10° psi and 40.83 x 10° psi. 
Aluminum and Aluminum only, becomes an acceptable solution. 
This hereby reconfirms the design already established in 
previous chapters that Aluminum 7075 should be used in 
Configuration F with a beam width of 1 inch. The units of 
M Hlasticity are *10° psi and the units for density are 


lb/in*, 


Contour Plot of Stiffness, Mass 


Density = 0.163397 
M_Elasticity = 14.9610 
Stiffness = 72004.4 
Mass = 11.4101 





Density = 0.101211 


M_Elasticity = 10.3442 
Stiffness = 48030.3 
| Mass = 7.05390 





Figure 32. MINITAB results of holding the beam width at 0.25 
inches. 
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Contour Plot of Stiffness, Mass 





Density = 0.163152 
M_Elasticity = 14.9510 
Stiffness 















Density = 0.100966 
M_Elasticity = 10.1635 







Density = 0.0654315 
M_Elasticity = 6.91136 
Stiffness = 70621,3 

Mass = 4.76231 


Figure 33. MINITAB results of holding the beam width 
at 0.5 inches. 


Contour Plot of Stiffness, Mass 





Figure 34. MINITAB results of holding the beam width 
at 1 inch. 
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VI. COMPOUND LOAD ANALYSIS 


A. ADDITIONAL ANALYSIS 


Now that the optimum material is determined to be 
space-grade Aluminum and the optimum configuration is the 
pattern designated by F, the next step can be taken. A 
thorough analysis needs to be done on Configuration F. The 
first analysis was conducted with the gravitational forces 
being applied in the negative Y direction. However, for this 
two-dimensional analysis, the structure also needs to be 
analyzed for forces in the positive and negative xX 
directions, the positive Y direction, and a combination of 


the two. Figures 35-38 show all force combinations analyzed. 


i It r 


600 Ib 600 Ib 600 Ib 600 Ib 600 Ib 600 Ib 
— > —_> <— <—- 
600 Ib OO lb 600 Ib 00 Ib 
—>*x —>X —>x 
F1 Fla Fib 
Figure 35. Diagrams of compound forces analyzed (F1I-F1b). 
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r r 


600 Ib 600 Ib 600 Ib 600 Ib 600 Ib 600 Ib 


| | f f I I 
600 Ib OO lb 600 Ib 00 Ib 
=x —>x —>4k 
Fic Fid Fle 


Figure 36. Diagrams of compound forces analyzed (Flc-Fle). 


600 Ib 600 Ib 600 Ib 
| , | 
600 Ib 600 Ib 
600 Ib 600 Ib 
—xX ——> —>x <— —>Xx 
Fif F1 Fih 
600 Ib ee 8 600 Ib 
Figure 37. Diagrams of compound forces analyzed (FIf-Flh). 


600 Ib 600 Ib 600 Ib 
I I 
600 Ib 600 Ib 
600 Ib 600 Ib 
—>x —-> —>xX <— —>X 
i j Fik 
coo b | Pa | mu etoi | 
Figure 38. Diagrams of compound forces analyzed (F11i-F1k). 


B. MATLAB RESULTS 


The compound loads from Figures 35-38 are put into 
MATLAB to return the maximum stress and deflections for each 
design. Appendix E contains the MATLAB code used for 
analysis, and Appendix F shows the stress and displacements 


for each element and node. The maximum values were then used 





in Appendix G to solve for how each combined load affects 


the reduction in stiffness. Figure 39 shows each setup's 





corresponding stress and deflection. From this chart, it can 
be seen that the largest deflection and stress is seen in 
setups a, b, d and e. Figure 40 shows how each setup 
compares in stiffness percentages. The ultimate stress of 
Aluminum 7075 is 68 ksi, and the yield stress is 57 ksi. The 
maximum setup stress is 6.2 ksi, which is vastly smaller 


than the limits of the material. 
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Figure 39. Stress and deflections for all setups. 
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Figure 40. Percent reduction in stiffness variations for 
compound loadings. 
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VII . CONCLUSION 


A. CONFIGURATION 


Mass is the leading cause of high cost prices for the 
construction and launch of satellites. Something with a lot 


of mass costs more in materials and takes more time to 





integrate all the systems. This integration also contributes 





a large portion to the full cost. As a result, the goal for 





any satellite is to have minimal mass. The NPSCuL structure 
carries up to 24 CubeSats into space. With the decrease in 
size of components for these satellites, the CubeSat 


designers are able to put more and more mass into the 10 cm? 





volume. These designers constantly desire more and more mass 
to be allowed. Optimizing the weight of NPSCuL can either 


allow more mass to be put into the CubeSats, or can go 





toward attempting to lower the magnification of forces seen 








at resonance frequencies during launch of the rocket that 





carries the structure to space. Looking at ways to optimize 
mass utilization of modifications to the modifiable area is 


ascertained. Looking at various truss like structures 








contained inside the MA to remove as much mass as possible, 
while maintaining some stiffness, was completed. Using an 


“X” shape of crossbeams in the desired area proved to be the 





best design. 


B. ALUMINUM 7075 AND COMBINED LOADING 


The NPSCuL structure does not have to operate in space 





for the designed life that a satellite does. It is merely 





the transportation of the satellites on their way to space. 
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As such, the next step in optimization was to evaluate 





whether Aluminum 7075 (space-grade Aluminum) was the best 
choice for materials. Analyzing the range of metals from 
Magnesium to Tungsten to Gold, a simultaneous analysis was 


conducted using the “X”-shaped configuration. The percent 





reduction in stiffness and the percent reduction in mass 
were evaluated at different crossbeam widths for all the 


metals using an analysis program called MINITAB. The 








constraint is that the structure cannot exceed 177 lbs. It 


was found that the crossbeam width must be 1 inch and the 





material must be an Aluminum alloy to meet the requirements. 





The initial analysis used only forces in one direction 





for calculations and comparisons. The next step to evaluate 





was if the compound loading, which the structure will see, 





renders the design a failure. The yield stress of Aluminum 
7075 is 57 ksi. A combination of directions and forces were 
applied and the maximum stress seen in any of the members 


was 6.2 ksi. 


Cc. FUTURE WORK 


This optimization of mass looked at the two-dimensional 
response of the separate walls. Each wall only had six ksi 
of a stress when analyzed with compound loads and the yield 
of the material is 57. The ultimate of the material is 68 


ksi and stresses seen wer ten times less than that. The 





next step in the process will be to take the modifications 





in the modifiable area and construct a model in IDEAS. From 


there, a dynamic response can be determined; if it is a 





plausible solution, the structure should be built and 








vibration tested. The analysis of nonmetal options for the 








structure should also be undertaken. 
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APPENDIX 


FEMA1 


nel=4; 

nnel=2; 

ndof=2; 

nnode=4; 
sdof=nnode*ndof; 


gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =4; 
nodes (3,1)=2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =4; 
bedof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bedof (3) =7; 

beval (3) =0; 

bedof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


££ (4) =-600; 
ff£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
Halt 


end 


[kk, 


nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











Tr 














k=FETRUSS2 (el, leng, area,0,beta,1); 


kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 




















k=FETRUSS2 (el, leng, area,0,beta,1); 


for i=1: (nnel*ndof) 
eldisp(i)=disp(index(i)); 
end 
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elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMA5 


nel=4; 
nnel=2; 
ndof=2; 
nnode=4; 


sdof=nnode* 


ndof; 


gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord (3,1)=6.365; gcoord (3,2) 
gcoord (4,1)=6.365;gcoord (4,2) 
elprop(1,1)=10300000; 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop(2,2)=.125; 

elprop (3,1)=10300000; 

elprop (3,2)=.125; 
elprop (4,1) =10300000; 

elprop (4,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =4; 
nodes (3,1)=2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =4; 
bcdof (1) =1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof, sdof) ; 
index=zeros (nnel*ndof,1); 


elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 


stress=zeros(nel,1); 


ff£ (4) =-600; 
f£f£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord (nd ( 


1),1);yl=gcoord(nd(1),2); 
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end 


[kk, 


x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 





Tr. 











k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1) *2+(y2-y1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 


55 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEMA25 


nel=4; 
nnel=2; 
ndof=2; 
nnode=4; 


sdof=nnode* 


ndof; 


gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord (3,1)=6.365;gcoord (3,2) 
gcoord (4,1)=6.365;gcoord (4,2) 
elprop (1,1)=10300000; 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 

elprop (2,2)=.0625; 

elprop (3,1)=10300000; 

elprop (3,2)=.0625; 
elprop (4,1) =10300000; 

elprop (4,2)=.0625; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =4; 
nodes (3,1)=2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =4; 
bcdof (1) =1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof, sdof) ; 
index=zeros (nnel*ndof,1); 


elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 


stress=zeros(nel,1); 


ff£ (4) =-600; 
f£f£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord (nd ( 


1),1);yl=gcoord(nd(1),2); 
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end 


[kk, 


x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 





Tr. 











k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1) *2+(y2-y1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 
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if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEMB1 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop(1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop (2,2)=.25; 

elprop (3,1)=10300000; 
elprop (3,2)=.25; 
elprop (4,1) =10300000; 
elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 
elprop (6,2)=.25; 

elprop (7,1)=10300000; 
elprop(7,2)=.25; 

elprop (8,1)=10300000; 
elprop (8,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =6; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2) =5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMB5 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop (1,1)=10300000; 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop (2,2)=.125; 

elprop (3,1)=10300000; 
elprop (3,2)=.125; 

elprop (4,1)=10300000; 
elprop (4,2)=.125; 

elprop (5,1)=10300000; 
elprop(5,2)=.125; 
elprop (6,1) =10300000; 
elprop (6,2)=.125; 

elprop (7,1)=10300000; 
elprop(7,2)=.125; 

elprop (8,1)=10300000; 
elprop (8,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =6; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2) =5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMB25 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop (1,1)=10300000; 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 
elprop(2,2)=.0625; 

elprop (3,1)=10300000; 
elprop (3,2)=.0625; 

elprop (4,1)=10300000; 
elprop (4,2)=.0625,; 

elprop (5,1)=10300000; 
elprop(5,2)=.0625; 

elprop (6,1)=10300000; 
elprop (6,2)=.0625,; 

elprop (7,1)=10300000; 
elprop(7,2)=.0625; 

elprop (8,1)=10300000; 
elprop (8,2)=.0625; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =6; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2) =5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMC1 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop(1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop (2,2)=.25; 

elprop (3,1)=10300000; 
elprop (3,2)=.25; 
elprop (4,1) =10300000; 
elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 
elprop (6,2)=.25; 

elprop (7,1)=10300000; 
elprop(7,2)=.25; 

elprop (8,1)=10300000; 
elprop (8,2)=.25; 

nodes (1,1)=1;nodes(1,2)=5; 
nodes (2,1)=1;nodes (2,2) =2; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2)=5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMC5 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop (1,1)=10300000; 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop (2,2)=.125; 

elprop (3,1)=10300000; 
elprop (3,2)=.125; 

elprop (4,1)=10300000; 
elprop (4,2)=.125; 

elprop (5,1)=10300000; 
elprop(5,2)=.125; 
elprop (6,1) =10300000; 
elprop (6,2)=.125; 

elprop (7,1)=10300000; 
elprop(7,2)=.125; 

elprop (8,1)=10300000; 
elprop (8,2)=.125; 

nodes (1,1)=1;nodes(1,2)=5; 
nodes (2,1)=1;nodes (2,2) =2; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2)=5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMC25 


nel=8; 

nnel=2; 

ndof=2; 

nnode=6; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =6.365;gcoord (4 
gcoord(5,1)=6.365;gcoord (5 
gcoord(6,1)=6.365; gcoord (6 
elprop (1,1)=10300000; 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 
elprop(2,2)=.0625; 

elprop (3,1)=10300000; 
elprop (3,2)=.0625; 

elprop (4,1)=10300000; 
elprop (4,2)=.0625,; 

elprop (5,1)=10300000; 
elprop(5,2)=.0625; 

elprop (6,1)=10300000; 
elprop (6,2)=.0625,; 

elprop (7,1)=10300000; 
elprop(7,2)=.0625; 

elprop (8,1)=10300000; 
elprop (8,2)=.0625; 

nodes (1,1)=1;nodes(1,2)=5; 
nodes (2,1)=1;nodes (2,2) =2; 
nodes (3,1)=2;nodes (3,2) =5; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2)=5; 
nodes (6,1) =3;nodes (6,2) =4; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
bcdof (1)=1; 

beval (1) =0; 

bcedof (2) =2; 

beval (2) =0; 

bedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 
kk=zeros (sdof,sdof); 


index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1)/(x2-x1)); 
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end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMD1 


nel=11; 

nnel=2; 

ndof=2; 

nnode=7; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) 
gcoord(5,1)=6.365;gcoord (5,2) 
gcoord(6,1)=6.365;gcoord (6,2) 
gcoord(7,1)=3.1825; gcoord(7,2 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 
elprop (2,1) =10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 

elprop (7,1)=10300000; 
elprop(7,2)=.25; 

elprop (8,1)=10300000; 

elprop (8,2)=.25; 

elprop (9,1)=10300000; 

elprop (9,2)=.25; 

elprop (10,1)=10300000; 
elprop(10,2)=.25; 

elprop (11,1)=10300000; 
elprop(11,2)=.25; 

nodes (1,1)=1l;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =7; 
nodes (3,1) =2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =7; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =7; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =7; 
nodes (9,1)=5;nodes (9,2) =6; 
nodes (10,1)=6; nodes (10,2) =7; 
nodes (11,1)=7;nodes (11,2)=1; 
bedof (1) = 
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ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1); y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











r 











k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 
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end 


num= 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMD5 


nel=11; 

nnel=2; 

ndof=2; 

nnode=7; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) 
gcoord(5,1)=6.365;gcoord (5,2) 
gcoord(6,1)=6.365;gcoord (6,2) 
gcoord(7,1)=3.1825; gcoord(7,2 
elprop(1,1)=10300000; 
elprop(1,2)=.125; 
elprop (2,1) =10300000; 
elprop(2,2)=.125; 

elprop (3,1)=10300000; 

elprop (3,2)=.125; 
elprop (4,1) =10300000; 

elprop (4,2)=.125; 

elprop (5,1)=10300000; 
elprop(5,2)=.125; 

elprop (6,1)=10300000; 

elprop (6,2)=.125; 

elprop (7,1)=10300000; 
elprop(7,2)=.125; 
elprop (8,1) =10300000; 

elprop (8,2)=.125; 

elprop (9,1)=10300000; 

elprop (9,2)=.125; 

elprop (10,1)=10300000; 
elprop(10,2)=.125; 

elprop (11,1)=10300000; 
elprop(11,2)=.125; 

nodes (1,1)=1l;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=7; 
nodes (3,1) =2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =7; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1) =4; nodes (6,2) =7; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =7; 
nodes (9,1)=5;nodes (9,2) =6; 
nodes (10,1)=6; nodes (10,2) =7; 
nodes (11,1)=7;nodes(11,2)=1; 
bedof (1) = 
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ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1); y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











r 











k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 
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end 


num= 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMD25 


nel=11; 

nnel=2; 

ndof=2; 

nnode=7; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) 
gcoord(5,1)=6.365;gcoord (5,2) 
gcoord(6,1)=6.365;gcoord (6,2) 
gcoord(7,1)=3.1825; gcoord(7,2 
elprop (1,1)=10300000; 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 

elprop (2,2)=.0625; 

elprop (3,1)=10300000; 

elprop (3,2)=.0625; 

elprop (4,1)=10300000; 

elprop (4,2)=.0625; 

elprop (5,1)=10300000; 

elprop (5,2)=.0625; 

elprop (6,1)=10300000; 

elprop (6,2)=.0625; 

elprop (7,1)=10300000; 
elprop(7,2)=.0625; 

elprop (8,1)=10300000; 

elprop (8,2)=.0625; 
elprop (9,1) =10300000; 

elprop (9,2)=.0625; 

elprop (10,1)=10300000; 

elprop (10,2)=.0625; 

elprop (11,1)=10300000; 
elprop(11,2)=.0625; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=7; 
nodes (3,1) =2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =7; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1) =4; nodes (6,2) =7; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =7; 
nodes (9,1)=5;nodes (9,2) =6; 
nodes (10,1)=6; nodes (10,2) =7; 
nodes (11,1)=7;nodes (11,2)=1; 
bedof (1) = 
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ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1); y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











r 











k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


85 


end 


num= 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEME1 


nel=12; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 





ae ON Ole ll 





gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord (3,2) 
gcoord (4,1)=3.1825; ae 
gcoord(5,1)=6.365; gcoord ( 
gcoord(6,1)=6.365; gcoord ( 
gcoord(7,1)=6.365;gcoord ( 
gcoord(8,1)=3.1825;gcoord 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 
elprop (2,1) =10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 
elprop (3,2)=.25; 
elprop (4,1) =10300000; 
elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 
elprop (6,1) =10300000; 
elprop (6,2)=.25; 

elprop (7,1)=10300000; 
elprop(7,2)=.25; 

elprop (8,1)=10300000; 
elprop (8,2)=.25; 

elprop (9,1)=10300000; 
elprop (9,2)=.25; 

elprop (10,1) =10300000 
elprop(10,2)=.25; 

elprop (11,1)=10300000; 
elprop(11,2)=.25; 

elprop (12,1)=10300000; 
elprop (12,2)=.25; 

nodes (1,1)=1;nodes(1,2)= 
nodes (2,1)=2;nodes(2,2)= 
nodes (3,1) =2;nodes (3,2)= 
nodes (4,1) =2;nodes (4,2)= 
nodes (5,1) =3;nodes(5,2)= 
nodes (6,1) =4;nodes (6,2)= 
nodes (7,1) =4;nodes(7,2)= 
nodes (8,1)=5;nodes (8,2)= 
nodes (9,1) =6;nodes (9,2)= 
nodes (10,1)=6;nodes (10,2 
nodes (11,1)=7;nodes (11,2 


nodes (12,1)=8; nodes (12,2)=1; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-400; 
ff£(8)=-400; 
££ (10) =-400; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 





leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








r 














k=FETRUSS2 (el, leng, area,0,beta,1); 


kk=FEASMBL1 (kk, k, index) ; 





end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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eae 


end 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-y1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEME5 


nel=12; 

nnel=2; 

ndof=2; 

nnode=8; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)= 
gcoord(2,1)=0;gcoord(2,2)= 
gcoord(3,1)=0;gcoord(3,2)= 
gcoord (4,1) =3.1825; gcoord ( 
gcoord(5,1)=6.365;gcoord(5, 
gcoord(6,1)=6.365;gcoord(6, 
gcoord(7,1)=6.365;gcoord (7 
gcoord (8,1)=3.1825; gcoord ( 
elprop (1,1)=10300000 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop(2,2)=.125; 

elprop (3,1)=10300000; 
elprop (3,2)=.125; 
elprop (4,1) =10300000; 
elprop (4,2)=.125; 

elprop (5,1)=10300000 
elprop(5,2)=.125; 

elprop (6,1)=10300000; 
elprop (6,2)=.125; 

elprop (7,1)=10300000 
elprop(7,2)=.125; 

elprop (8,1)=10300000; 
elprop (8,2)=.125; 
elprop (9,1) =10300000; 
elprop(9,2)=.125; 

elprop (10,1) =10300000 
elprop(10,2)=.125; 

elprop (11,1)=10300000; 
elprop (11,2)=.125; 

elprop (12,1)=10300000; 
elprop(12,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =8; 
nodes (3,1)=2;nodes (3,2) =4; 
nodes (4,1) =2;nodes (4,2) =3; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1) =4; nodes (6,2) =6; 
nodes (7,1)=4;nodes (7,2) =5; 
nodes (8, 1)=5;nodes (8,2) =6; 
nodes (9, 1)=6;nodes (9,2) =8; 
nodes (10,1)=6;nodes (10,2)= 
nodes (11,1)=7;nodes (11,2)= 


O; 

6; 

Ls 
4,2)=12; 
2)=12; 
2)=6; 
72) =07 
8,2) =6; 
ee 

8; 


nodes (12,1)=8; nodes (12,2)=1; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-400; 
ff£(8)=-400; 
££ (10) =-400; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 





leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) / (x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








r 














k=FETRUSS2 (el, leng, area,0,beta,1); 


kk=FEASMBL1 (kk, k, index) ; 





end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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eae 


end 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-y1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEME25 


nel=12; 
nnel=2; 
ndof=2; 
nnode=8; 
sdof=nnode*n 


gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 





elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 11, 1 
Li7:2 
12/4 
12,2 





l= 
) 
= 
) 
ME 
) 


dof; 


=0;gcoord(1, 
=0;gcoord(2, 
=0;gcoord (3,2 
-1825;gcoord 
- 365; gcoord ( 
- 365; gcoord ( 
- 365; gcoord ( 
-1825;gcoord 


2)5 
2) = 
) 


a NOUR | 


10300000; 
=.0625; 

10300000; 
=.0625; 

10300000; 
=.0625; 


O; 

6; 

Ls 
4,2)=12; 
,2)=12; 
,2)=6; 
72) =07 
8,2) =6; 
ee 

8; 


nodes (12,1)=8; nodes (12,2)=1; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


ff£(6)=-400; 
ff£(8)=-400; 
££ (10) =-400; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 





leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








r 














k=FETRUSS2 (el, leng, area,0,beta,1); 


kk=FEASMBL1 (kk, k, index) ; 





end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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eae 


end 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-y1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF1 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (4) =-600; 
ff£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 











kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 


97 





Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF'5 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord (5,2) =6; 
elprop (1,1)=10300000; 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop(2,2)=.125; 

elprop (3,1)=10300000; 

elprop (3,2)=.125; 
elprop (4,1) =10300000; 

elprop (4,2)=.125; 

elprop (5,1)=10300000; 
elprop(5,2)=.125; 

elprop (6,1)=10300000; 
elprop(6,2)=.125; 

elprop (7,1)=10300000; 
elprop(7,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bcdof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (4) =-600; 
ff£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 











kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
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Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF25 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 

elprop (2,2)=.0625; 

elprop (3,1)=10300000; 

elprop (3,2)=.0625; 

elprop (4,1)=10300000; 

elprop (4,2)=.0625; 

elprop (5,1)=10300000; 

elprop (5,2)=.0625,; 

elprop (6,1)=10300000; 

elprop (6,2)=.0625; 

elprop (7,1)=10300000; 

elprop (7,2)=.0625,; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (4) =-600; 
ff£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 











kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
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Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMG1 


nel=13; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 


gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 





elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 


cop 
cop 
cop 
cop 





=0;gcoord(1,2)= 
=0;gcoord(2,2)= 
=0;gcoord (3,2) 
- 365; gcoord ( 
. 365; gcoord ( 
- 365; gcoord ( 
-1825;gcoord 
-1825;gcoord 


AA OUS I 


=10300000; 
S2o5F 
=10300000; 
s297 
=10300000; 
E297 
=10300000; 
w2oF 
=10300000; 
250% 
=10300000; 
Zor 
=10300000; 
sD Dy 
=10300000; 
3207 
=10300000; 


=10300000; 
10300000; 
~ 25; 


Io We th Wek thet 
rn 


omowo es CO W OO ~~) N 


Ne 
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10,1)= 
Lad; 
12,1 


L3¢d 


10,2 
ely ad 
L252 
13,2 


;nodes 


nodes 5 
=5;nodes 
6 
7 


nodes 
nodes 
nodes 


, 


;nodes 
;nodes 


( ) 7 

( ) =6; 
( ) 7; 
( ) 1 


aN 


) 
) 
) 
= ) 


a 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 














k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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tor 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEMG5 


nel=13; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) =12; 
gcoord(5,1)=6.365;gcoord (5,2) =6; 
gcoord(6,1)=6.365;gcoord (6,2) =0; 
gcoord(7,1)=3.1825; gcoord(7,2)=3; 
gcoord(8,1)=3.1825; gcoord (8,2) =9; 
elprop (1,1)=10300000 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 

elprop (2,2)=.125; 

elprop (3,1)=10300000; 

elprop (3,2)=.125; 
elprop (4,1) =10300000; 

elprop (4,2)=.125; 

elprop (5,1)=10300000 
elprop(5,2)=.125; 

elprop (6,1)=10300000; 

elprop (6,2)=.125; 

elprop (7,1) =10300000 
elprop(7,2)=.125; 

elprop (8,1)=10300000; 

elprop (8,2)=.125; 
elprop (9,1) =10300000; 
elprop(9,2)=.125; 

elprop (10,1)=10300000; 
elprop(10,2)=.125; 

elprop (11,1)=10300000; 
elprop(11,2)=.125; 

elprop (12,1)=10300000; 

elprop (12,2)=.125; 

elprop (13,1)=10300000; 

elprop (13,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 

nodes (2,1)=2;nodes (2,2) =7; 

nodes (3,1)=2;nodes (3,2) =8; 

nodes (4,1) =2;nodes (4,2) =3; 

nodes (5,1)=3;nodes (5,2) =8; 

nodes (6,1) =3;nodes (6,2) =4; 

nodes (7,1)=4; nodes (7,2) =8; 

nodes (8, 1)=4; nodes (8,2) =5; 

nodes (9,1)=5;nodes (9,2) =8; 
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10,1)= 
Lad; 
12,1 


L3¢d 


10,2 
ely ad 
L252 
13,2 


;nodes 


nodes 5 
=5;nodes 
6 
7 


nodes 
nodes 
nodes 


, 


;nodes 
;nodes 


( ) 7 

( ) =6; 
( ) 7; 
( ) 1 


aN 


) 
) 
) 
= ) 


a 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 














k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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tor 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEMG25 


nel=13; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) =12; 
gcoord(5,1)=6.365;gcoord (5,2) =6; 
gcoord(6,1)=6.365;gcoord (6,2) =0; 
gcoord(7,1)=3.1825; gcoord(7,2)=3; 
gcoord(8,1)=3.1825; gcoord (8,2) =9; 
elprop (1,1)=10300000 

elprop (1,2)=.0625,; 

elprop (2,1) =10300000 
elprop(2,2)=.0625; 

elprop (3,1) =10300000 

elprop (3,2)=.0625,; 

elprop (4,1) =10300000 

elprop (4,2)=.0625; 

elprop (5,1)=10300000 
elprop(5,2)=.0625; 

elprop (6,1)=10300000 

elprop (6,2)=.0625,; 

elprop (7,1) =10300000 
elprop(7,2)=.0625; 

elprop (8,1) =10300000 

elprop (8,2)=.0625,; 

elprop (9,1) =10300000 

elprop (9,2)=.0625,; 

elprop (10,1)=10300000; 

elprop (10,2)=.0625; 

elprop (11,1)=10300000; 
elprop(11,2)=.0625; 

elprop (12,1) =10300000; 

elprop (12,2)=.0625; 

elprop (13,1)=10300000; 

elprop (13,2)=.0625; 

nodes (1,1)=1;nodes (1,2) =2; 

nodes (2,1)=2;nodes (2,2) =7; 

nodes (3,1)=2;nodes (3,2) =8; 

nodes (4,1) =2;nodes (4,2) =3; 

nodes (5,1)=3;nodes (5,2) =8; 

nodes (6,1) =3;nodes (6,2) =4; 

nodes (7,1)=4; nodes (7,2) =8; 

nodes (8, 1)=4; nodes (8,2) =5; 

nodes (9,1)=5;nodes (9,2) =8; 
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10,1)= 
Lad; 
12,1 


L3¢d 


10,2 
ely ad 
L252 
13,2 


;nodes 


nodes 5 
=5;nodes 
6 
7 


nodes 
nodes 
nodes 


, 


;nodes 
;nodes 


( ) 7 

( ) =6; 
( ) 7; 
( ) 1 


aN 


) 
) 
) 
= ) 


a 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
ff£(8)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 














k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
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tor 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:l1:nel; 
stresses=[numm' stress] 
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FEMH1 


nel=14; 
nnel=2; 
ndof=2; 
nnode=8; 
sdof=nnode*n 


gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 
gcoor 





elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 
elpr 


cop 
cop 
cop 
cop 
cop 
cop 





dof; 


=0;gcoord(1,2)= 
=0;gcoord(2,2)= 
=0;gcoord (3,2) 
-365;gcoord 
-365;gcoord 
.365;gcoord 
.1825;gcoor 
-1825;gcoord 


( 
( 
( 
d 


AA OUS I 


=10300000; 


S2o5F 


=10300000; 


s297 


=10300000; 


E297 


=10300000; 


w2oF 


=10300000; 


250% 


=10300000; 


Zor 


=10300000; 


sD Dy 


=10300000; 


3207 


=10300000; 


=10300000; 
=10300000; 
10300000; 
2297 
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nodes (8, 1)=4; nodes (8,2) =8; 
nodes (9,1)=4;nodes (9,2) =5; 
nodes (10,1)=5; nodes (10,2) =8; 
nodes (11,1)=5;nodes (11,2) =7; 
nodes (12,1)=5;nodes (12,2) =6; 
nodes (13,1)=6; nodes (13,2) =7; 
nodes (14,1)=7;nodes (14,2)=1; 
bedof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
f£(8)=-600; 
for iel=l:nel 


end 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord ( 
x2=gcoord ( 


leng=sqrt ( ( 


nd(1),1);yl=gcoord(nd(1),2 
nd(2),1);y2=gcoord(nd(2),2 


*%2-x1)*2+(y2-yl1)%2); 


if (x2-x1)==0; 


beta=2* 
else 


atan(1); 


beta=atan ((y2-yl1) /(x2-x1)); 


end 


el=elprop(iel,1); 
area=elprop(iel,2); 











[index]=FE 


ELDOF (nd,nnel,ndof); 





k=FETRUSS2 ( 





kk=FEASMBL1 








1,leng,area,0,beta,1); 


(kk, k, index) ; 
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[kk, 


ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


ee 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMH5 


nel=14; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) =12; 
gcoord(5,1)=6.365;gcoord (5,2) =6; 
gcoord(6,1)=6.365;gcoord (6,2) =0; 
gcoord(7,1)=3.1825; gcoord(7,2)=3; 
gcoord(8,1)=3.1825; gcoord (8,2) =9; 
elprop (1,1)=10300000 
elprop(1,2)=.125; 

elprop (2,1)=10300000; 
elprop(2,2)=.125; 

elprop (3,1)=10300000; 

elprop (3,2)=.125; 
elprop (4,1) =10300000; 

elprop (4,2)=.125; 

elprop (5,1)=10300000 
elprop(5,2)=.125; 

elprop (6,1)=10300000; 

elprop (6,2)=.125; 

elprop (7,1) =10300000 
elprop(7,2)=.125; 

elprop (8,1)=10300000; 

elprop (8,2)=.125; 
elprop (9,1) =10300000; 
elprop(9,2)=.125; 

elprop (10,1)=10300000; 
elprop(10,2)=.125; 

elprop (11,1)=10300000; 
elprop(11,2)=.125; 

elprop (12,1) =10300000; 
elprop(12,2)=.125; 

elprop (13,1)=10300000; 

elprop (13,2)=.125; 

elprop (14,1) =10300000; 

elprop (14,2)=.125; 

nodes (1,1)=1;nodes (1,2) =2; 

nodes (2,1)=2;nodes (2,2) =7; 

nodes (3,1)=2;nodes (3,2) =5; 

nodes (4,1) =2;nodes (4,2) =8; 

nodes (5,1)=2;nodes (5,2) =3; 

nodes (6,1) =3;nodes (6,2) =8; 

nodes (7,1) =3;nodes (7,2) =4; 
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nodes (8, 1)=4; nodes (8,2) =8; 
nodes (9,1)=4;nodes (9,2) =5; 
nodes (10,1)=5; nodes (10,2) =8; 
nodes (11,1)=5;nodes (11,2) =7; 
nodes (12,1)=5;nodes (12,2) =6; 
nodes (13,1)=6; nodes (13,2) =7; 
nodes (14,1)=7;nodes (14,2)=1; 
bedof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
f£(8)=-600; 
for iel=l:nel 


end 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord ( 
x2=gcoord ( 


leng=sqrt ( ( 


nd(1),1);yl=gcoord(nd(1),2 
nd(2),1);y2=gcoord(nd(2),2 


*%2-x1)*2+(y2-yl1)%2); 


if (x2-x1)==0; 


beta=2* 
else 


atan(1); 


beta=atan ((y2-yl1) /(x2-x1)); 


end 


el=elprop(iel,1); 
area=elprop(iel,2); 











[index]=FE 


ELDOF (nd,nnel,ndof); 





k=FETRUSS2 ( 





kk=FEASMBL1 








1,leng,area,0,beta,1); 


(kk, k, index) ; 
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[kk, 


ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


ee 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMH25 


nel=14; 
nnel=2; 
ndof=2; 
nnode=8; 


sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0; gcoord (2,2) =6; 
gcoord(3,1)=0;gcoord(3,2)=12; 
gcoord (4,1)=6.365;gcoord (4,2) =12; 
gcoord(5,1)=6.365;gcoord (5,2) =6; 
gcoord(6,1)=6.365;gcoord (6,2) =0; 
gcoord(7,1)=3.1825; gcoord(7,2)=3; 
gcoord(8,1)=3.1825; gcoord (8,2) =9; 
elprop (1,1)=10300000 
elprop(1,2)=.0625; 

elprop (2,1)=10300000; 

elprop (2,2)=.0625; 
elprop (3,1) =10300000; 

elprop (3,2)=.0625; 

elprop (4,1) =10300000 

elprop (4,2)=.0625; 

elprop (5,1)=10300000 

elprop (5,2)=.0625; 

elprop (6,1)=10300000 

elprop (6,2)=.0625; 

elprop (7,1) =10300000 
elprop(7,2)=.0625; 

elprop (8,1) =10300000 

elprop (8,2)=.0625; 

elprop (9,1) =10300000 

elprop (9,2)=.0625; 

elprop (10,1)=10300000; 
elprop(10,2)=.0625; 

elprop (11,1)=10300000; 
elprop(11,2)=.0625; 

elprop (12,1) =10300000; 

elprop (12,2)=.0625; 

elprop (13,1)=10300000; 

elprop (13,2)=.0625; 

elprop (14,1) =10300000; 

elprop (14,2)=.0625; 

nodes (1,1)=1;nodes (1,2) =2; 

nodes (2,1)=2;nodes (2,2)=7; 

nodes (3,1) =2;nodes (3,2) =5; 

nodes (4,1) =2;nodes (4,2) =8; 

nodes (5,1)=2;nodes (5,2) =3; 

nodes (6,1) =3;nodes (6,2) =8; 

nodes (7,1)=3;nodes (7,2) =4; 
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nodes (8, 1)=4; nodes (8,2) =8; 
nodes (9,1)=4;nodes (9,2) =5; 
nodes (10,1)=5; nodes (10,2) =8; 
nodes (11,1)=5;nodes (11,2) =7; 
nodes (12,1)=5;nodes (12,2) =6; 
nodes (13,1)=6; nodes (13,2) =7; 
nodes (14,1)=7;nodes (14,2)=1; 
bedof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcedof (3)=11; 

beval (3) =0; 

bcdof (4) =12; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof); 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


f£(6)=-600; 
f£(8)=-600; 
for iel=l:nel 


end 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord ( 
x2=gcoord ( 


leng=sqrt ( ( 


nd(1),1);yl=gcoord(nd(1),2 
nd(2),1);y2=gcoord(nd(2),2 


*%2-x1)*2+(y2-yl1)%2); 


if (x2-x1)==0; 


beta=2* 
else 


atan(1); 


beta=atan ((y2-yl1) /(x2-x1)); 


end 


el=elprop(iel,1); 
area=elprop(iel,2); 











[index]=FE 


ELDOF (nd,nnel,ndof); 





k=FETRUSS2 ( 





kk=FEASMBL1 








1,leng,area,0,beta,1); 


(kk, k, index) ; 
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[kk, 


ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


ee 


end 


num= 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 





stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 


if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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APPENDIX B 


7.0000 0003 = 7.0000 G@oook)§= 70000 goood )=s 7.0000 0.0000) 7.0000 acoD = 7.0000 0.0000 
80000 -0.0046 80000 -0.0023 BOND -0.0112 00ND -0.0056 8.0n00 -0.0028 40000 -A.0086 
9.0000 0.0006 9.0000 0.0008 9.0000 0.0053 9.0000 0.0026 9.0000 0.0013 9.0000 0.0000 


10.0000 -0.0023 «610.0000 -0.0011 100000 0006 100000 02028 100000 02004 10000 

10000 -0.0000 8110000 0.0000 110000 G0000) 6110000 0.0000 110000 @.0000 §=61L0000 0.00065 
120000 0.0000) =— 120000 ao T2.0000 -00000 §=120000 0.0000 120000 -20000 12.0000 -0.005 
13.0000 -0.0000 8130000 0.0000 130000 O0000) 86130000 00000 13.0000 a0000) = 13.0000 0000 
14,0000 -0.0008 «= 14.0000 -0.0004 14.0000 -00000 14.0000 0.0000 14.0000 -00000 14.0000 0.0000 
15.0000 3.0000 


15.0000 0 15.0000 0 15.0000 0 
16.0000 -0.0036 16.0000 -0.0018 16.0000 -0.0112 16.0000 -0.0056 16,0000 -0.0028 16.0000 a.com 


streses= stresses = stresses = stresses= stresses = streses= 

1L0e+08 * L0e+008 * L0e-+008 * 1L.0e+008 * L0e+008 * L.0e+003 * 

Qn 38798 Como 1999 CONN 9.6000 oni 4.8000 aomin 2.4000 agmio 9.6000 oom 4.8000 gam 2.4000 
ano -125 foo -16eq7 onon amm = Ao oom = 6oo~e anm = oam~ 1.673 Onn ase 

0.0030 18782 0.0080 09141 00030 -0.0000 0.0030 -2.0000 0.0080 0.0000 0.0030 -3.4223 


0.0040 -11709 = 0.0080 -05855 0040 9.6000 = @.0080 4.8000 n0n% 24000 0.00% 6.4000 
0.0050 3.9968 005) 1.9981 aD 
0.0060 -11709 0.0060 -05855 0008 00000 
00070 08571 OOo 0420 aon”) -0.0000 
0.0080 -11709 0.0080 -0.5855 0080 9.6000 


0.0090 3.9968 0.0080 1.9984 aon o0000 0.0080 0.0090 ©0000 0.0090 16973 

0.0100 -1175 00100 -05855  O0100 -00000 O0100 -20000 O0100 0.0000 anole 9.6000 ono 43000 
O.0110 -13415 O0110 0607 a0110 9.6000 ao10 4.8000 O00 24000 ao10 00000 0.0110 -c.0mD 
00120 388 86001 1999 O01 00ND §8=6 i 0D OD) ODS. ONO 
0.0130 -13415 0.0130 0607 00130 -00000 00130 20000 0.0130 -0.0000 

Q.0140 -13415 0.0140 -0.677 


1.0000 =O iO i1o000 Clio iLomo0 600 1.0000 «=O 10000 «sO iLoodo )6=Olio 

12.0000 «OO 129000:—=—oO 1000: =O 12.0000 :~«~O 120000:«O 120000 «OO 
stresses = stresses = stresses = stresses = stresses = stresses = 
1.0e+003 * 1.0e1003 * L.0e+003 * L.0e+003* 1L.0e1003 * ‘L0e+003 * 


ae ee a a ae 
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APPENDIX C 


Fetaterfal properties ‘rom [19] 


F (gfem“8} €(Gpal leeen'e 

Pall cd eel lM ame 
p.50-090 0 ft .ssfo.2es) | peoatolanse] hoosen[se72s) pt 
Se ae a 
ec pet bea bff J 
FeaironGredesoaois ——~i80800 Sb rooase) | 
enn eee fee Espa fa — at — fn. — fs 
Pluminum Aley Sse,cunsementTs [Lor _}#7_a.sefoosm} p2alios) _ps3__| 


173) 


1.77 [0.0889) | 


albania =f es — 
ade L 


Geld{fecldwerked) (4 800-10,280 [14100 lig.sztoee7] p72) paz os(so) _—bzorsay = 
Pemerestwened mora bastard bras] bse [besa] _ba 
Bilver(cokdwerked) ft7O2ia00 BY io.g9os79) patio.) sr | Sees tas) = 
jae fom Be pov besa bat hemi bate lo 
Menalgog SONGS B29 be.sof0.218) leofze) psa _—eo(ss) _—ssoyeay to 
TaD ATOR] — evaas — fia — eee peel — pe — pet eel 
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APPENDIX D 
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APPENDIX 


FEMF la 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 





gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord (3,1)=6.365; gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord (5,2) =6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 
elprop (2,1) =10300000; 

elprop (2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 

elprop (7,1)=10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2) =5; 
nodes (3,1) =2;nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bcdof (2) =2; 

beval (2) =0; 

bedof (3) =7; 

beval (3) =0; 

bedof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof,sdof) ; 
index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
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eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 


££ (3) =600; 
££ (4) =-600; 
££ (5) =600; 
ff£(6)=-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord (nd (1) 
x2=gcoord(nd(2),1); y2=gcoord (nd (2) 


leng=sqrt ( (x2-x1)%*2+(y2-y1)%*2); 


if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1)/(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 








Tr 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 





[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 
disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1); yl=gcoord (nd (1) 
x2=gcoord (nd(2),1); y2=gcoord (nd (2) 


leng=sqrt ( (x2-x1)*2+(y2-y1)%*2); 


if (x2-x1)==0; 
beta=2*atan(1); 
else 
beta=atan ((y2-yl1) /(x2-x1)); 
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p2)ye 
r2)3 


r2)7 
p2)7 


end 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd, nnel,ndof) ; 














k=FETRUSS2 (el, leng, area,0,beta,1); 








for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2)%*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF 1b 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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ff ( 
fete s( 
ff ( 
EE i 


for 


end 


[kk, 


)=-600; 
)=-600; 
)=-600; 
)=-600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF 1c 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (4) =600; 
££ (6) =600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 











kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
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Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF1d 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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ff ( 
ie ie( 
ff ( 
EE ( 


for 


end 


[kk, 


)=600; 
)=600; 
)=600; 
)=600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


141 


end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 


142 


FEMF le 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=3;nodes (5,2) =4; 
nodes (6,1)=4; nodes (6,2) =5; 
nodes (7,1)=5;nodes(7,2)=1; 
bcdof (1)=1; 

beval (1) =0; 

bedof (2) =2; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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ff ( 
fete ( 
ff ( 
EE ( 


for 


end 


[kk, 


)=-600; 
)=600; 
)=-600; 
)=600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF1f 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2) =1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (2)=-600; 
££ (4) =-600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 











kk=FEASMBL1 (kk, k, index) ; 


end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
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Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF1g 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2)=1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ ( 
Ee '( 
ff ( 
ff ( 


for 


end 


[kk, 


)=600; 
)=-600; 
)=600; 
)=-600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF 1h 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2)=1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ ( 
ee "( 
ff ( 
ff ( 


for 


end 


[kk, 


)=-600; 
)=-600; 
)=-600; 
)=-600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF 11 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2)=1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ (2)=600; 
ff (4) =600; 


for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 
beta=2*atan(1); 

else 
beta=atan ((y2-yl1)/(x2-x1)); 

end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 











= 








k=FETRUSS2 (el, leng, area,0,beta,1); 








kk=FEASMBL1 (kk, k, index) ; 





end 


[kk, ££] =FEAPLYC2 (kk, ££, bcdof,bcval); 





disp=kk\ff; 
for iel=l:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
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Tr 





‘-ELDOF (nd, nnel,ndof) ; 





index=E 








k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp(i)=disp (index (i)); 
end 





elforce=k*eldisp; 
stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce(3))<0; 
stress (iel)=-stress(iel); 
end 





end 


num=1:1:sdof; 
displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF15 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2) =1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ ( 
Ee ( 
ff ( 
ff ( 


for 


end 


[kk, 


1 
Zz 
3 
4 


)=600; 
)=600; 
)=600; 
)=600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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FEMF 1k 


nel=7; 

nnel=2; 

ndof=2; 

nnode=5; 
sdof=nnode*ndof; 








gcoord(1,1)=0;gcoord(1,2)=0; 
gcoord(2,1)=0;gcoord(2,2)=12; 
gcoord(3,1)=6.365;gcoord (3,2) =12; 
gcoord (4,1)=6.365;gcoord (4,2) =0; 
gcoord(5,1)=3.1825; gcoord(5,2)=6; 
elprop (1,1)=10300000; 
elprop(1,2)=.25; 

elprop (2,1)=10300000; 
elprop(2,2)=.25; 

elprop (3,1)=10300000; 

elprop (3,2)=.25; 
elprop (4,1) =10300000; 

elprop (4,2)=.25; 

elprop (5,1)=10300000; 
elprop(5,2)=.25; 

elprop (6,1)=10300000; 

elprop (6,2)=.25; 
elprop (7,1) =10300000; 
elprop(7,2)=.25; 

nodes (1,1)=1;nodes (1,2) =2; 
nodes (2,1)=2;nodes (2,2)=5; 
nodes (3,1) =2; nodes (3,2) =3; 
nodes (4,1) =3;nodes (4,2) =5; 
nodes (5,1)=4;nodes (5,2) =5; 
nodes (6,1)=4;nodes (6,2) =1; 
nodes (7,1)=5;nodes (7,2)=1; 
bcdof (1)=5; 

beval (1) =0; 

bcdof (2) =6; 

beval (2) =0; 

bcdof (3) =7; 

beval (3) =0; 

bcdof (4) =8; 

beval (4) =0; 


ff=zeros(sdof,1); 

kk=zeros (sdof, sdof) ; 

index=zeros (nnel*ndof,1); 
elforce=zeros (nnel*ndof,1); 
eldisp=zeros (nnel*ndof,nnel*ndof) ; 
stress=zeros(nel,1); 
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££ ( 
Ee ( 
ff ( 
ff ( 


for 


end 


[kk, 


)=-600; 
)=600; 
)=-600; 
)=600; 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+ (y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 


el=elprop(iel,1); 
area=elprop(iel,2); 
[index]=FEELDOF (nd, nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





kk=FEASMBL1 (kk, k, index) ; 





ff£]=FEAPLYC2 (kk, f£,bcdof,bcval) ; 





disp=kk\ff; 


for 


iel=1:nel 


nd(1)=nodes(iel,1); 
nd(2)=nodes(iel,2); 


xl=gcoord(nd(1),1);yl=gcoord(nd(1),2); 
x2=gcoord(nd(2),1);y2=gcoord(nd(2),2); 


leng=sqrt ((x2-x1) *2+(y2-yl1)%2); 





if (x2-x1)==0; 

beta=2*atan(1); 
else 

beta=atan ((y2-yl1) /(x2-x1)); 
end 
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end 


num= 


el=elprop(iel,1); 
area=elprop(iel,2); 
index=FEELDOF (nd,nnel,ndof) ; 

















k=FETRUSS2 (el, leng, area,0,beta,1); 





for i=1: (nnel*ndof) 
eldisp (i)=disp (index (i)); 
end 





elforce=k*eldisp; 


stress (iel)=sqrt (elforce (1) *2+elforce (2) %*2) /area; 





if ((x2-x1) *elforce (3) )<0; 
stress (iel)=-stress(iel); 
end 





1:1:sdof; 


displ=[num' disp] 


numm=1:1:nel; 
stresses=[numm' stress] 
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APPENDIX F 


Stress and deflections from compound loading of Configuration F 


FL 


stresses = 


1.0e+003 * 


0.0010 114485") 0.0010 


10.0020 
10.0030 
10.0040 
10.0050 
10.0060 
10.0070 


Flh 
displ = 


1.0000 
2.0000 
3.0000 
4.0000 
5.0000 
‘6.0000 
7.0000 
8.0000 
9.0000 


1.0759 
0.5041 
-1.0759 


1.0759 
1.0759 


0.0006 
-0.0043 
-0.0022 


0.0000 
-0.0000 
-0.0000 
-0.0000 
-0.0005 


10.0000 -0.0020 


stresses = 


1.0e+003 * 


0.0010 
0.0020 
10.0030 
10.0040 
10.0050 
0.0060 
10.0070 


0.3216 
-2.3527 


3.0807 
-2.3527 
0.9564 
3.0807 


Fla 


stresses = 


1.0e+003 * 


0.0020 
0.0030 
0.0040 
0.0050 
0.0060 
0.0070 


Fli 


displ = 


1.0000 
2.0000 
3.0000 
4.0000 
5.0000 
6.0000 
7.0000 
8.0000 


3.0752 
0.5041 
4.0459 
5.9743 


4.0459 


1.2730 


Fib 


1.0000 0.0006 
2.0000 0.0043 
3.0000 0.0022 


5.0000 -0.0000 
6.0000 0.0000 
7.0000 0.0000 
8.0000 0.0000 
9.0000 0.0005 
10.0000 0.0020 


stresses = 


1.0e+003 * 


0.0010 0.3216 
0.0020 2.3527 
0.0030 

0.0040 -3.0807 
0.0050 2.3527 
0.0060 -0.9564 


27167) 0.0070 -3.0807 


Fic 


stresses = 


1,0e+003 * 


stresses = 


1.0e+003 * 


0.0010 4485") 0.0010 5.9743 


0.0020 1.0759 
0.0030 -0.5041 
0.0040 1.0759 
0.0050 

0.0060 1.0759 
0.0070 1.0759 


Flk 


displ = 


1.0000 -0.0022 


0.0020 -4.0459 
0.0030 -0.5041 
0.0040 
0.0050 3.0752 
0.0060 -4.0459 
0.0070 


2.0000 0.0047 


3.0000 -0.0006 
4.0000 0.0043 
5.0000 0.0000 
6.0000 -0.0000 
7.0000 0.0000 
8.0000 0.0000 
9.0000 -0.0005 
10.0000 0.0020 


stresses = 


1.0e+003 * 


0.0010 0.3216 
G.0020 3.0807 
9.0030 0.9564 
0.0040 -2.3527 
0.0050 3.0807 
0.0060 

0.0070 -2.3527 
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1.0000 -0.0000 
2.0000 -0.0000 
3.0000 -0.0242 
4.0000 -0.0036 


6.0000 0.0070 
7.0000 0.0000 
8.0000 -0.0000 
9.0000 -0.0072 
10.0000 0.0008 


stresses = 


1.0e+003 * 


0.0010 3.0752 
0.0020 

0.0030 -0.5041 
0.0040 -4,0459 
0.0050 5.9743 
0.0060 

0.0070 -4.0459 


Fif 
displ = 


1.0000 0.0008 
2.0000 

3.0000 -0.0008 
4.0000 -0.0085 
5.0000 0.0000 
6.0000 -0.0000 
7.0000 0.0000 
8.0000 0.0000 
9.0000 -0.0000 
10.0000 -0.0020 


stresses = 
1.0e+003 * 


0.0010 «(0 
0.0020 

0.0030 1.2730 
0.0040 

0.0050 

0.0060 -1.2730 
0.0070 


Fig 
displ = 


1.0000 0.0022 
2.0000 

3.0000 0.0006 
4.0000 -0.0043 
5.0000 -0.0000 
6.0000 0.0000 
7.0000 -0.0000 
8.0000 -0.0000 
9.0000 0.0005 
10.0000 -0.0020 


stresses = 
1.0e+003 * 


0.0010 0.3216 
0.0020 -3.0807 
0.0030 -0.9564 
0.0040 2.3527 
0.0050 -3.0807 
0.0060 

0.0070 2.3527 
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APPENDIX G 







































































F=0*A 
Absolute | Absolute | g (psi) 
Max Max A k=F/5 

deflection Stress (in’2) (Ibs) (Ibs/in) 
F1 0.0017 1.4495 | 1449.5 0.25 | 362.375 | 213,162 
Flal 0.0245 6.1978 | 6197.8 0.25 | 1549.45 63,243 
F1b1 0.0245 6.1978 | 6197.8 0.25 | 1549.45 63,243 
Fic1 0.0017 1.4495 | 1449.5 0.25 | 362.375 | 213,162 
F1d1 6197.8 0.25 |} 1549.45 63,243 
Fle1 6197.8 0.25 | 1549.45 63,243 
F1f1 0.0045 2.7167 | 2716.7 0.25 | 679.175 | 150,928 
Flg1 0.0047 3.5024 | 3502.4 0.25 875.6 | 186,298 
Fih1 0.0047 3.5024 | 3502.4 0.25 875.6 | 186,298 
F1i1 0.0045 2.7167 | 2716.7 0.25 | 679.175 | 150,928 
F1j1 0.0047 3.5024 | 3502.4 0.25 875.6 | 186,298 
875.6 | 186,298 
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Pcr (Ibs) Status 
919.2069 Safe 
919.2069 
919.2069 
919.2069 Safe 
919.2069 
919.2069 
919.2069 Safe 
919.2069 Safe 
919.2069 Safe 
919.2069 Safe 
919.2069 Safe 
919.2069 Safe 
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